Shape resonances in modified effective-range theory for electron-molecule collisions 
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We develop a simple model of shape resonances in electron-molecule collisions that is based on the 
modified effective-range expansion and analytical solutions of the Schrodinger equation for the long- 
range part of the interaction potential. We apply our model to electron scattering on N2 and CO2. 
The parameters of the effective-range expansion (i.e. the scattering length and the effective range) 
are determined from experimental, integral elastic cross sections in the 0.1 - 1.0 eV energy range. 
For both molecular targets our treatment predicts shape resonances that appear slightly higher than 
experimentally known resonances in total cross sections. Agreement with the experiment can be 
improved by assuming the position of the resonance in a given partial wave. Influence of quadrupole 
potential on resonances is also discussed: it can be disregarded for N2 but gets significant for CO2. 
In conclusion, our model developed within the effective range formalism reproduces well both the 
very low-energy behavior of the integral cross section as well as the presence of resonances in the 
few eV range. 

PACS numbers: 34.80.Bm 
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I. INTRODUCTION 

2 II temporary negative ion states in N2 and CO2 at 
2.1 eV and 3.8 eV 0, Q respectively, are the best known 
examples of so-called shape resonances, in which the in- 
coming electron is captured within the potential well of 
the target, rather than to a specific electronic orbital. 
This distinction came clearly from the first discovery of 
the shape (in N2) and Fcshbach (in He y] ^resonances by 
Schulz, see detailed discussion in Ref. [J, [fa]. 

The widths and amplitudes of these resonances, as ap- 
pearing in total cross sections, are well established, see 
for ex. Ref. 0,0,9 for N 2 and [UQ3 for C0 2 . In both 
resonances the vibrational excitation contributes signif- 
icantly to the total cross section (about 1/6 in N 2 and 
1/3 in CO2, see [TH); both targets are non-polar, with 
negligible direct rotational excitation cross sections. 

However, apart from the resonances, the low-energy 
dependences of total cross section in N 2 and C0 2 dif- 
fer completely: the cross section rises sharply in the 
zero-energy limit in C0 2 [12] and falls slowly in N 2 Q. 
Particularly, the C0 2 scattering received a rich theo- 
retical treatment, starting from the early assignment of 
this state as shape 2 n„ state [l~3[ and calculations us- 
ing model potentials, for ex. by Morrison and co-workers 
14, [l5|, More recent works, both experimental [l7| 
and theoretical [l8[ , indicate that the mechanism of this 
resonance is complex, with coupling between different vi- 
brational modes; the analysis requires the vibrational ex- 
citations to be included [19j. However, with few excep- 
tions [HI, H(| calculations were performed separately for 
the zero-energy and resonance energy ranges, see for ex. 

In our previous work [25| we revisited a rather old prob- 
lem - application of the modified effective range theory 



(MERT) to electron and positron scattering [2fJ . In that 
theory the scattering is due to the polarization potential 
and the short range potential is modelled by an effective 
range expansion. These features make MERT particu- 
larly suitable for the near-to-zero en ergy range, treated 
with difficulties in other approaches |27l [28j . The ncar- 
to-zero energy MERT expansion continues to be applied 
to numerous problems: integral cross sections for atoms 
[29l . l30j , non-polar [3l|, [32| and polar molecules [H, H|j] . 
MERT is frequently used by experimentalists for extrapo- 
lation of measured differential and integral cross sections 
towards energies and angles inaccessible by experiment 
HI, H, H3- Different extensions of MERT, consisting 
in developing phase shifts (or elements of the reactance, 
K matrix) into series of k were applied by Morrison et 
al. [27| to rotational excitation of N 2 molecule, and by 
Macri and Barrachina to electron - metastable He scat- 
tering and electron photodetachment from Li - [38l . |39| . 

The most extensive, to our knowledge, analysis of 
MERT applicability in electron scattering was done for 
noble gases by Buckman and Mitroy [401 ] . They con- 
cluded that MERT fails above about 1 eV, exact energy 
depending on the target. However, in that as well in 
numerous other MERT applications, those were integral 
cross sections which were developed directly in series of 
the colliding electron momentum k. As we showed in the 
previous work [25j . the range of applicability of such a 
series is limited to a few tenths of eV. Therein [25| we 
proposed a novel way of applying MERT. It consists in 
calculating the phase shifts due to the polarization poten- 
tial from Mathieus solutions of the Schrodinger equation 
and introducing the effective-range expansion exlusively 
for the short-range potential. There is no explicit k ex- 
pansion for the integral cross sections in our approach; 
only the short-range contribution expressed in terms of 
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an appropriate boundary condition, is expressed as se- 
ries of k. This makes the numerical procedures more 
complicated but allows to extend the applicability of 
MERT for positron-argon and -nitrogen scattering to the 
range of 2-3 eV. We note that our method closely resem- 
bles quantum-defect-theory formulation of the scatter- 
ing problem, where the quantum-defect parameter, that 
slowly varies with energy, is expanded in a series of k 
(see e.g. [H|). In electron-atom scattering the quantum- 
defect theory has been already applied by Watanabe and 
Greene [42j to analyze photodetachment process of K~ 
and by Fabrikant [431 ] to low-energy electron-metal atoms 
resonances. 

In this paper we apply our MERT-bascd model to the 
two molecules N2 and CO2, well studied by beam and 
swarm techniques at energies down to less than 0.1 eV. 
The parameters of MERT expansion for the short-range 
potential are derived from fitting procedure using the rec- 
ommended [44] integral cross sections below 1 cV. Next, 
we use these potentials to calculate cross sections up to 
10 eV. In both molecules broad p-wavc resonances ap- 
pear. In such an unconstrained fit the resonances appear 
at higher energies than experimental values. By small, 
within the experimental error bar, modifications of the 
data used for the inversion procedure, one can reproduce 
the exact position of the resonances. In that case both 
resonances appearing in the p or d-partial wave channels 
are narrower than experimentally determined (but this 
can be partially caused by neglecting the nuclear motion, 
see for ex. Ref. [l!| ) . We discuss also the relative contri- 
butions to integral cross sections from the non-spherical 
part of the polarizability and the quadrupole moment of 
molecule, using the distorted-wave approximation. 

The paper is organized as follows. In section[TT]we dis- 
cuss the analytical solutions of the Schrodinger equation 
with the polarization potentials, and we present princi- 
ples of the modified effective-range expansion. In sec- 
tion [TTT] we demonstrate that MERT accounts for shape 
resonances in p and d waves, assuming some simplified 
model of electron-molecule potential. The actual electron 
resonances in N2 and CO2 are analyzed in section IIVI 
Section [V] discusses correction to MERT cross sections 
due to the non-spherical part of the polarization poten- 
tial and the molecule quadrupole moment. We present 
some conclusions in section [VTl Finally, three appendices 
give some technical details of the analytical solutions for 
the polarization potential. 



II. SCATTERING OF A CHARGED PARTICLE 
ON A MOLECULE AND THE 
QUANTUM-DEFECT APPROACH 

The relative motion of a charged particle and a non- 
polar molecule is described by the Schrodinger equation 



where p, is the reduced mass, E is the relative energy of 
the particles, and V(v) is the particle- molecule potential. 
V(r) can be written as a sum 



V(r) = V 1 (r) + V 2 {r) + V s (r) 
of long-range contributions 
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(3) 
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and a short-range potential Vs(r). The long-range po- 
tential is expressed in terms of the static spherical and 
non-spherical polarizabilities a and a 2 , respectively, and 
the quadrupole moment Q. The short-range part Vs(r) 
comes into play at the distances comparable to the size 
of the molecule, where the molecule cannot be treated 
as a single object. Numerous approaches were proposed 
in the past, based on modifications of the polarization 
potential [45|], deriving the short-range interaction from 
electronic densities etc. [461. o r expanding the scattering 
amplitude at low energies [431 ] . In the present approach 
we express Vs(r) in terms of some boundary conditions 
imposed on the wave function at r — > 0, while the effects 
resulting from the finite range of V s (r) are included ex- 
plicitly in the framework of the modified effective-range 
theory. 

In the following we assume that the nonisotropic part 
V^r) can be neglected in comparison to the isotropic 
polarization potential Vi(r). This is justified as long as 
the quadrupole moment Q and the nonsphcrical polariz- 
ability a 2 expressed in atomic units are much smaller 
than ao. In such a case the main effect of the non- 
isotropic terms is the coupling between different partial 
waves, while their contribution to the total cross-section 
at low energies remains small. As we show later in sec- 
tion Sec. [Vj where we analyze the corrections arising 
from V 2 (r) , this is a good approximation for N2 and CO2 
molecules. The MERT analysis for these molecules is 
performed in Sec. IIVI 

The radial part of the Schrodinger equation with the 
isotropic polarization potential V\ (r) reads 



d 2 _ 1(1 + 1) (R*) 2 2[iE 
dr 2 r 2 r 4 K 2 



*!(»•) = 0, (5) 



- — A + V(r)-E 



*(r) = 0, 



(1) 



where vE'i(r) denotes the radial wave function for the par- 
tial wave /. For the polarization potential, it is con- 
venient to introduce some characteristic units R* and 
E* , where R* = y/ ae 2 /j,/h 2 denotes the characteristic 
length of the r~ 4 potential, and E* = h 2 /(2fiR* 2 ) is the 
characteristic energy. With appropriate change of vari- 
ables, the Schrodinger equation §5§ can be transformed 
into Mathieu's differential equation of the imaginary ar- 
gument [1^, [47], EH , and solved analytically in terms of 
the continued fractions. Some basic properties of the an- 
alytical solutions are discussed in Appendices El [B] and 



[Cl (see [49|, [5(| for more details on Mathieu functions). 
Here, we focus only on the behavior at small and large 
distances. For r <C R*, when the polarization potential 
dominates over centrifugal potential and the constant en- 
ergy term, behavior of *&i(r) is given by 



/' R* 

*i( r ) rsin ( h (!>i 

1 r 



(6) 



where <pi is a short-range phase, which is determined by 
V s (r). For E = and / = the solution ^ becomes 
exact at all distances, and from its asymptotic behavior 
at large distances one can easily determine the value of 
the s-wave scattering length 



-R* cot( 



(7) 



At large distances: r>i!*, ^i{r) must take the form of 
the scattered wave 



^i(r) r ~°° sin(fcr -l%+r)l), 



(8) 



where k = y/2fiE/h. Using the Mathieu functions one 
can find the following relation between the phase shift r\\ 
and the short range phase <pi [26[ 



tan rji 



m 2 - tanS 2 + tan(</>; + Z§) tan<5(m 2 - 1) 
tan 6(1 - m 2 ) + tan(0 ; + Z§)(1 - m 2 tan 2 6) ' 

(9) 

where 5 = ^(u — Z — i), and m and v are some parame- 
ters, that are determined by the analytical solutions (see 
Appendices lAl and IBI for details). 

In general, the parameter <pi entering the asymptotic 
formula ©, depends on energy, and can be expanded in 
powers of k. In our case it is more convenient to expand 
tan(0; + Zf ), entering formula ©: 



tan(0 ; + Zf ) = A t + -R*Rik 2 



(10) 



where Ai = tan(0; +Z§)| - The lowest order correc- 
tion in k is quadratic, and can be in terp reted as an effec- 
tive range Ri for the partial wave Z [26j . 

Combining ©, (flO|l . (|Cljl . (|C2]> one can easily obtain 
the low-energy expansion of the phase shifts [2y| 
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(12) 



Here, g = fci?*, a = a s /R* = -1/A , and 6 = A\. The 
low-energy expansions (p~X|) - (fT2|) are applicable only at 
small energies J? « e*. In contrast, in our approach we 
use the exact formula © for the phase shift, with the 
short-range phase expanded according to (jTT 
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FIG. 1: The model potential Eq. (|13p that is used to test 
MERT for the shape resonances 



III. MODEL CALCULATION 

To illustrate how MERT describes shape resonances for 
the long-range r -4 potential we consider a simple model 
of the potential with a square-well interaction at short 
distances (see FigfT]) 



V(r) 



-U 



r < R m 
R m <r<R 
r > Ro 



(13) 



Here, R m denotes a hard-core radius and Ro is the radius 
of the square-well potential. The potential V(r) does not 
reproduce the full behavior of the electron- molecule inter- 
action, however, it captures the basic properties of shape 
resonances in such a system, and allows to verify the ac- 
curacy of our MERT-based approach. The Schrodingcr 
equation for the potential V(r) and partial-wave quan- 
tum number Z can be solved exactly: 



'0, r<R m 
*(r) = { rji(xr) + Cryi(xr), 
Ay/TM u (ln 



•Jkr 



R m <r<R a 
B V^(ln^ 



r>Ro 
(14) 

where H 2 x 2 / (2/i) = E + U and constants A, B, C are to 
be determined from the continuity conditions for W(r) at 
r = R m and r = Rq. With the help of the asymptotic 
formulas (|B3f and (|E4]) for M v (z) and T u (z) for large 
and negative arguments (large r) we find the phase shifts 
and the total cross-section for scattering on V(r). We 
compare the exact results with our MERT model where 
the MERT coefficients Ai and Ri are found by expanding 
the formula tan(0/ + v\ + |) = A/B (c.f. Eq. (|B6]l) in 
the powers of g, and utilizing series expansion (|C1[) for 
the characteristic exponent. 

Fig. [2] shows the dependence of the total elastic cross- 
section on energy, for few particular values of U, for which 
p-wave shape-resonances appear. The range Rq was cho- 
sen to be smaller than R*, but still, of the same order 
as R* , that roughly corresponds to the conditions of the 
low-energy electron(positron)-molecule scattering. The 
resonances show up as peaks at some particular values 
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FIG. 2: Total elastic cross-section versus energy for the scat- 
tering in the model potential (|13[) , calculated for R m — O.li?*, 
Ro = 0.5R* and for different depths U of the square-well po- 
tential. Results are scaled by the the characteristic distance 
R* and the characteristic energy E* of the polarization poten- 
tial. The peaks correspond to the p-wave shape resonances, 
that are accurately reproduced by MERT. 



FIG. 3: Same as Fig. f2] but for values of U for which shape 
resonances in d wave appear. 



IV. MERT ANALYSIS OF LOW-ENERGY 
ELECTRON SCATTERING 



of the energy. We observe that MERT, including only 
two lowest expansion coefficients A[ and Ri of the short- 
range potential, reproduces very accurately the exact re- 
sults. The accuracy decreases for higher energies, where 
one can observe some discrepancies, in particular small 
peaks corresponding to the g-wave resonances that ap- 
pear at different energies for exact and MERT curves. 
This results from higher sensitivity of the quasi-bound 
states on the parameters of the short-range potential for 
larger /, and the absence of the higher order terms in 
expansion (|10j) when E becomes too large. Nevertheless, 
the approximate MERT treatment remains very accurate 
for the main resonances in p-wave channel, even for en- 
ergies significantly larger than E* . 



A similar comparison, but for <i-wave resonances is 
presented in Fig. [3] fn this case the resonance peaks 
are narrower due to the stronger centrifugal barrier, 
and weaker coupling between quasi-bound and scatter- 
ing states. Again, the agreement between MERT and 
the exact results is very good, and some discrepancies 
can be observed only at E ~ 10E*. Finally, we note 
that for I = 2 and parameters of Fig. [3] the quasi-bound 
state is totally localized within the range of the short- 
range potential. This shows that MERT is able to repro- 
duce shape resonances resulting solely from the short- 
range potential, even if the latter is treated in terms of 
an effective-range expansion. 



We investigate the low-energy electron scattering on 
N2 and CO2 molecules using a theoretical model based 
on MERT for s and p partial waves [25[ . Our model con- 
tains four unknown parameters: the scattering length 
a = a s /R* = — 1/A and the effective range i? for s 
wave, and the zero-energy contribution Ai and the effec- 
tive range R\ for p wave, that are determined by fitting 
to the experimental data. In this way for I = and 
/ = 1 we neglect corrections of the order higher than q 2 
in expansion ()10|) . On the other hand, for partial waves 
with Z > 2 we include only the lowest-order contribution 
to the scattering phase shifts that are due to the long- 
range polarization potential (fT2"| . This is sufficient as 
long as the energy E remains smaller than the centrifu- 
gal barrier, that for I = 2 has height of 9E* . The total 
elastic cross-section a{k) = Airk~ 2 J ^ ll sin 2 rji(k) is deter- 
mined from the phase shifts © and (fT2"|) . with v and m 
evaluated numerically, using the procedure described in 
Appendices O and [B] 

In the calculations we use recent experimental values 
of the polarizability, as measured in electron scattering 
experiments: a = 11.54ao 3 for N2 and a = 16.92ao 3 for 



CO2 |5f|. First, we fit our model by the least-square, not 



weighted method to the experimental data exclusively 
below the resonances, i.e. for energies E < 1.2eV and 
E < 2eV for N2 and CO2 , respectively, down to the low- 
est energies available from beam experiments [44T |. The 
fitted parameters from this check are listed in Table [T] 
Note, that the scattering potential parameters are ex- 
pressed in characteristic distance units. For example, 
in N2 the scattering length a s amounts to 0.404 atomic 
units, close to the value of 0.420 a.u. derived by Morri- 
son et al. [27|. Subsequently, we extend MERT analysis, 



5 




-i-i 1 1 — ■ — i — i — i—i-i-j 1 1 — i — i — i — i i 1 1 r i i i i 

0,1 1 10 0,1 1 

Electron energy (eV) Electron energy (eV) 



FIG. 4: Total cross-section for the scattering of electrons 
on N2 versus the energy. Depicted are: recommended experi- 
mental data from review [iil ] (stars) , the theoretical fits based 
on MERT, assuming the resonance at 2.1eV in s wave (dot- 
ted line), p wave (dashed line), and d wave (dot-dashed line), 
and without assumption with respect to the position of the 
resonance (solid line) . The inset shows in addition the s-wave 
and p-wave contributions to the MERT fit not assuming the 
position of resonance. 





R*(a ) 


E*(eV) 


a 


Ai 


Ro/R* 


Ri/R* 


N 2 


3.397 


1.179 


0.119 


-0.537 


0.110 


0.262 


C0 2 


4.113 


0.8041 


-1.606 


-3.537 


-0.325 


1.208 



TABLE I: Characteristic distance R* , characteristic energy 
E* and four fitting parameters: a = a s /R* (s-wave scattering 
length), A\ (zero-energy contribution from the short-range 
potentials for p-wave), Ro (s-wave effective range), and Ri 
(p-wave effective range) for N2 and CO2 in the case of uncon- 
strained fits and using only the very low energy experimental 
data. 



using the fitted parameters to higher energies. Figs. 2] 
and [5] compare the experimental data (points) with the 
theoretical curves obtained from this fit (solid lines). For 
both molecules resonance maxima appear in the scatter- 
ing cross-section, however they are located at somewhat 
higher energies that the experimental ones (about 6 eV 
for N2 and 5 eV for CO2). For both molecules these 
"spontaneously" appearing resonances are due to the p- 
wave, see inserts in Fig. 2] and El In the case of CO2 the 
p-wave is subject to rapid changes also below the reso- 
nance. This is to be attributed to large values of b and 
R\ of the potential (see table UJ). Note that in experiment 
[52| some vibrational modes (001) show resonant maxima 
below the main 3.6 eV peak. 

The range of fitted data corresponds to quite low ener- 
gies in terms of characteristic units: E < E* (E < 2.5E*) 
for N2 (C02). Therefore, the fitting procedure, within 
the error bar of experimental points, is not able to pre- 



FIG. 5: The same as fig. [4] but for the scattering of electrons 
on CO2 with the recommended experimental data from review 
[■441 ]. Note, missing curve corresponding to s-wave resonance, 
which does not fit well to the experimental data in this case. 



diet with sufficient accuracy the cross section in the en- 
ergy range corresponding to resonance peaks. Also, ne- 
glection of the nonisotropic part of the potential leads 
to some additional errors. Moreover, our scheme ne- 
glects the possibility of d-wave or higher partial wave 
resonances, that can be suggested e.g. by the behav- 
ior of differential cross-section for the scattering on N2. 
To eliminate these shortcomings we have subsequently 
added to the fitting procedure also some experimental 
points, located in the energy range above the resonance: 
6eV < E < lOeV (N 2 ), and 5.5eV < E < 7eV (C0 2 ). In 
addition we have considered three possible scenario as- 
suming that resonances occur in partial waves s, p or d, 
exactly at the energy of the resonance peak observed in 
the experimental data. This corresponds to the require- 
ment that the phase shift r\i assumes value tt/2 at the 
peak location for I = 0,1,2, respectively. In this way we 
obtain some constrains on the six fitting parameters: a, 
A\, A\ (zero-energy contributions) and i?o, R\ and R2 
(effective ranges). The results are shown in Figs. [4] and 
presenting apart from the previous data, also theoretical 
curves with resonances assumed in s, p and d waves. The 
fitted parameters are listed in tables HT1 and IIIII 

We observe that for N2 the height of the resonance 
peak for s wave agrees well with the experimental data 
(we show experimental values averaged over the vibra- 
tional structure seen in the resonance, following recom- 
mended data from ref. |44|). This would suggest a possi- 
ble s-wave origin of the resonance. But this point requires 
some more detailed discussion. 

As said before, the generally accepted classification 
of resonances describes "shape" resonances as occurring 
within the shape of the effective potential, formed from 
the static and polarization contributions with the cen- 
trifugal barrier added. Obviously, for the s-wave the cen- 
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a Ax A 2 Ro/R" Ri/R* R 2 /R* 

s-wave res. -0.447 -0.021 - -3.100 -0.104 

p-wave res. 0.157 -4.314 -0.130 -47.05 4.704 0.079 

d-wave res. 0.123 -0.390 -0.222 -0.301 -0.174 0.293 



ferent molecular geometries, resulting from the bent con- 
figuration in the dominating vibrationally excited mode 
(010). 



TABLE II: Six fitting parameters: a = a s /R* (s-wave scat- 
tering length), A\ (zero-energy contribution for p-wave), A 2 
(zero-energy contribution for d-wave), Rq (s-wave effective 
range), Ri (p-wave effective range) and R2 (d-wave effective 
range) for N2 in the case of fits assuming shape resonances in 
s, p and d waves. The fits are made to the data below and 
above the resonance (see text for details). Note that for s- 
wave resonance only four parameters (two partial waves) were 
sufficient to fit the data. 





a 


Ai 


A 2 


Ro/R* 


Ri/R* 


R 2 /R* 


p-wave res. 


-1.616 


-6.250 


-0.406 


-0.266 


2.780 


0.131 


d-wave res. 


-1.605 


-2.470 


-2.624 


-0.340 


-0.234 


1.180 



TABLE III: Six fitting parameters: a — a s /R* (s-wave scat- 
tering length), A\ (zero-energy contribution for p-wave), A 2 
(zero-energy contribution for d-wave), Ro (s-wave effective 
range), Ri (p-wave effective range) and R 2 (d-wave effective 
range) for CO2 in the case of fits assuming shape resonances 
in p and d waves. The fits are made to the data below and 
above the resonance (see text for details). 



V. CORRECTIONS DUE TO THE 
NONISOTROPIC PART OF THE LONG-RANGE 
POTENTIAL 

Already early works stressed importance of other, 
apart from the sperical polarization, components of the 
long-range potential. For both molecules the non-sperical 
part of the polarizability is rather strong (see table IV); 
for CO2 also the quadrupole moment is significant. In 
the present analysis we approximate the quadrupole and 
non-spherical polarizability contributions together, but 
using two distinct methods for the zero-energy limit and 
at low E ~ E* energies. To analyze the contribution 
from the nonisotropic part of the potential at very low 
energies (E <C E*), we can compare the small-fc series 
expansions for the whole potential V\ (r) + V2 (r) [5(| , to 
the expansion including the spherical part V\ (r) only. For 
the isotropic part, the total cross-section for scattering on 
N2 and CO2 molecules can be obtained from expansions 
of the phase shifts given by Eqs. (llip - (fT2|) . Substituting 
the parameters of the potential listed in the Table HVl one 
obtains 



trifugal barrier does not exist so trapping of the electron 
within the shape of the effective potential is impossible. 
However, our calculations show a quick change in the 
s-wave phase shift, occurring approximately in the re- 
gion of the experimentally observed resonance so follow- 
ing the interpretation by Fano [H, |]34| we would classify 
this state as a resonance. On the other hand, assumption 
of resonances in p and d waves would lead to much nar- 
rower and higher peaks in the total cross section than the 
experimental values. Therefore, a simple interpretation 
of the the N2 resonance is not easy: i) clastic differen- 
tial cross sections at the resonance maximum (2.47 eV) 
show a dominating p-wave contribution, ii) differential 
cross sections for the first vibrational excitation (z/=0- 
1) at 2.45 eV show almost a perfect d-wave dependence 
[52j while iii) the MERT analysis does not exclude the 
resonance in the s-partial wave. Clearly, several partial 
waves contribute to the resonance and MERT analysis is 
not able to separate these contributions. 

For CO2 the MERT fit with constraints does not im- 
prove much the agreement with the experimental total 
cross sections, sec Fig. 3. The MERT analysis excludes 
completely the s-wave resonance, and both p and d wave 
resonances differ much in amplitude and width from the 
experimental values. For CO2 resonance the interplay of 
different partial waves is even more complex than for N2. 
On the other hand, recent calculations [H, H3| and ex- 
periments showed that the broad peak in CO2 cross 
section is rather a superposition of two resonances for dif- 



cn 2 {q) 



-in 



2A.60qa + 31.32g 2 \og(q)a 2 + 0(q 3 ) 



(15) 



&co 2 (q) 



4tt 



35.44ga + 45.12 9 2 log(q)a 2 + 0(q 3 ) 



(16) 
(17) 



On the other hand, application of the small-A; expansion 
for the full potential [56[ leads to 



°N 2 (q) =a 2 + 2 4.93 ga + o.23<7 + 0.11 



4?r 



Q"co 2 (g) 

47T 



+ 31.63(a - 0.15)(a + 0.16)g 2 log(g) + 0(q 3 ) 

(18) 

: a 2 + 39.56<?a- 3.97q + 1.32 

+ 49.09(a - 0.53)(a + 0.56)g 2 log(g) + 0(q 3 ) 

(19) 



We observe that for the typical scales of the scatter- 
ing length a ~ 1, the differences in the coefficients of 
Eas.([Tg ]) -(fT9"] ) and Eqs.|JTSJ)-|[TBJ arc at most of the order 
of 10%. 

For higher energies E ~ E* the importance of the 
nonisotropic part can be assessed with the help of the 
distorted wave approximation (57| . The scattering am- 
plitude can be written as a sum of the contributions from 
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a(fflp) a 2 ( a o) Q(eap) 



N 2 


11.54 


3.08 


-1.09 


C0 2 


16.92 


9.20 


-3.86 



TABLE IV: Parameters of the long-range potential VL(r) for 
molecules considered in our analysis. For Qo we use recent 
experimental values of [U, and for a 2 and Q we use exper- 
imental values used in Refs. [23| for N 2 and in Ref. [l4| for 

co 2 



E/E* erg R = 0.5R* Ro = 0.75JT R = R* 



1.0 


3.145 


0.0306 


0.0311 


0.0300 


2.0 


4.014 


0.0147 


0.0135 


0.0138 


4.0 


7.428 


0.0059 


0.0052 


0.0053 



TABLE V: The total elastic cross-section ao and corrections 
due to the nonspherical part of the long-range potential in 
N 2 molecule, for different energies E of the scattered electron 
and different short-range cut-off parameters Ro- 



V\ and y 2 separately (see e.g. [58|) 

/(k f ,kO = -^jL «kF|Vi|*+(ki)) + ($-(k f )|F 2 |* + )) , 

(20) 

where $ ± (k) is the scattering state for the isotropic part 
V 1 

$ ± (k) = e' kr + G±y 1 $ ± (k), (21) 
G^ is the free propagator 

Gq = Um(B - p 2 /2m± ie) -1 , (22) 

and is the scattering state for the full potential, sat- 
isfying the the following Lippmann-Schwinger equation 

= $+(ki)+G+V^ + , (23) 
Gt = \im(E - p 2 /2m — VI + ie)~ l . (24) 

The second term in Eq. (|20p can be written in the form 
of a Born series 

(*-(k f )|7 2 |* + ) = ($-(k f )|y 2 + V.G+V, + . . . |$ + (k;)). 

(25) 

Finally, the total cross-section can be calculated applying 
the optical theorem a = ^Im/(ki, k). 

After averaging over different orientations of molecules, 
the lowest order contribution from the nonisotropic part 
of the potential (<& _ (kf)|V2|$ + (ki)) vanishes, and for 
small anisotropic part, the leading contribution is pro- 
vided by ($-(k f )|V2G i "y 2 |* + (ki)). We evaluate this 
term using analytical solutions for the spherical poten- 
tial V*i and introducing an additional cut-off parameter 
i?0i that describes the distance when the interaction po- 
tential takes a form of the long range potential V\ +V2, 
which is roughly of the order of the size of the molecule. 

Tables [V] and IVII list the the corrections due to the 
nonisotropic part V2, together with total cross sections 
coming from the isotropic part V\. The analytical un- 
perturbed solution for V\ was calculated for parameters 
of Table (P), corresponding to the low-energy fits. For 
N2 the corrections are small, of the order of 1%. On 
the other hand, for CO2, the corrections are as large as 
40% of the total cross-section ao, therefore in this case 
an accurate analysis would require full inclusion of Vi- 
We note that corrections slightly depend on the cut-off 



E/E* go -Rp = 0.5R* Ro = 0.75R* R = R* 



1.0 


2.185 


0.174 


0.796 


0.698 


2.0 


1.345 


0.455 


0.514 


0.363 


4.0 


0.693 


0.207 


0.213 


0.153 



TABLE VI: The total elastic cross-section ao and corrections 
due to the nonspherical part of the long-range potential in 
C0 2 molecule, for different energies E of the scattered elec- 
tron and different short-range cut-off parameters Ro. 

parameter i?o, that is a consequence of the fact that the 
analytical solutions for the purely polarization potential 
V\ are not physical at small distances, where the short- 
range part of the potential comes into play. Therefore 
corrections listed in tables [V] and IVII estimate contribu- 
tions of V2 from large distances only, where the potential 
takes the asymptotic form V\ + V2 ■ 



VI. CONCLUSIONS 

MERT analytical solution was applied previously [25[ 
for positron scattering cross sections. Those cross sec- 
tions fall monotonically with energy, up to the positro- 
nium formation threshold, and do not exhibit particular 
structures, see (59j . The model approximated well Ar and 
N 2 integral elastic cross sections at 0-2 eV. The derived 
scattering potentials were characterized by negative val- 
ues of the scattering length-like parameters: a and A±, 
for s and p waves, respectively. 

At present, the MERT analysis was applied for electron 
scattering on non-polar molecules, N 2 and C0 2 . Param- 
eters of the short range potential (the scattering length 
and the effective range) for the s, p, and for higher ener- 
gies also the d partial wave, were fitted. The model ap- 
proximates very well the experimental data up to about 
1-2 eV. The s-wave scattering length deduced is posi- 
tive for N 2 and negative for C0 2 . Deriving the potential 
parameters from the experimental data below resonance 
and using them again for our model at higher energies, we 
obtained p-wave resonances for both molecules, although 
at slightly higher energies than observed experimentally. 
To check for the consistency of the fit, the model was 
applied jointly to the data below and above resonance, 
assuming that the resonance occur the in s, p or d wave, 
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at the energy as observed experimentally. Within the 
latter constraint, for N 2 a "resonance- like" peak, corre- 
sponding to the tt/2 phase change appears also in the 
s-wave channel. In that case the peak is not as sharp as 
pure shape resonances in p and d partial waves. 

As far as amplitude are considered, in N2 fixing the 
resonance position reproduces well the absolute value of 
the total cross section in its maximum, but only if the res- 
onance is assumed to appear in the s-wave channel. Such 
a fixed fit does not improve the agreement in the case of 
CO2. More experimental points at low energies, also for 
N 2 , would be useful for such a modelling By the way, 
we recall that Ramanan and Freeman [60j], on the basis 
of their swarm experiment even stated: "experimental 
data presently available do not rule out the existence of 
a Ramsauer-Townsend minimum at an electron energy of 
0.4 meV or lower" . 

Summarizing, the present work unifies in a single 
model the very low-energy dependence of the integral 
cross sections in N 2 and C0 2 with the occurrence of 
shape resonances at a few eV energies. It confirms that 
limitations on applicability of MERT to the very low en- 
ergy range, come rather from direct expansions of the 
phase-shifts in scries of fc, rather than from any intrinsic 
defects of the effective range model. Obviously, several 
aspects of resonant scattering, like partitioning between 
the elastic and vibrational/rotational excitation are be- 
yond the capacities of our model. It would be challenging 
to incorporate MERT potentials to more advanced the- 
oretical treatments, in particular of inelastic direct and 
resonant scattering (l^.l20Ll2l|. 



following form 0, E| 

OO 

M v (z) = (-l) n c n (^)J 2 „ + ,(2Vgcoshz),(A2) 

n— — 00 

00 

T„{z) = (-l)"c„(^)r 2 „ + ,(2Vgcoshz),(A3) 

n— — oo 

Here, v denotes the characteristic exponent, and J u (z) 
and Y v (z) are Bessel and Neumann functions respec- 
tively. Substituting the ansatz (fA2|) and ([A3]) into (jAlj) 
one obtains the following recurrence relation: 

[(2n + vf - a] c n + g(c n _! + c n+1 ) = 0, (A4) 

which can be solved in terms of continued fractions. 
To this end we introduce = c„/c„_i and h~ — 
c_„/c_„ + i for n > 0, which substituted into (|A4ft gives 
the continued fractions 



with d n = (2n + v) 2 — a. Characteristic exponent has 
to be determined from Eq. (|A4j) for n = expressed in 
terms of and hf . In practice, to find numerical values 
of the coefficients c n we set = and ft." = for some, 
sufficiently large m and calculate /i+ and h~ up to n= 1 . 
using (|A5j) . 

APPENDIX B: ASYMPTOTIC EXPANSIONS 
FOR LARGE ARGUMENTS 
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APPENDIX A: SOLUTIONS OF THE 
SCHRODINGER EQUATION FOR 1/r 4 
POTENTIAL 

The method for solving the Schrodingcr equation with 
the polarization potential was discussed in our previous 
paper (25[. Below we indicate the basic steps of the 
derivation. In the radial Schrodingcr equation ([5]) we 
substitute r = \f~R* e~ z j \fk and \&z(r) = ip{r)yJr/R*, 
which yields the Mathieu's modified differential equation 

mm 

d 2 ib 

— 4- - [a - 2q cosh 2z] i> = 0. (Al) 

dz l 

where a = (I + i) 2 and q = kR* . Two linearly indepen- 
dent solutions M(z) and T(z) can be expressed in the 



Asymptotic behavior of M v (z) and T v (z) for large z fol- 
lows immediately from asymptotic expansions of Bessel 
functions for large arguments [50| 

A^(z)^y|^ S! ,cos(e*V?-i^-i) (Bl) 

T v {z) Z =^\jl e -^^ sin (c*V5 - |^ - f ) (B2) 

where s v = X)^L-oo °n{ v )- To obtain asymptotic behav- 
ior for large and negative z it is necessary to connect 
solutions M v (z) and T v (z), with another pair of solution 
M v {— z) and T„(—z) across z = [48j |. In this way one 
obtains [25[ 

M v {z) z ^°^I^l mSl/ cos (^qe- z + \v-\) (B3) 

TM - ^ [sin (^qe-* + f v - f ) 

- cot ttv{itl 2 - 1) cos {^/qer z + \v - | ) 

(B4) 

where m = lim z ^ + M v {z)/M_ v (z). To evaluate m nu- 
merically, we rather avoid using formula (|A2[) . where 
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summation over n converges very slowly for z —> 0. An 
alternative approach is to use a different representation 
for the solutions of the Mathieu's equation J49| 



large distance described by Eqs. ([6])- (J9j) , can be readily 
obtained from asymptotic expansions (|Blj) - (|B4j) . 



W v {z) 



E 



(u)e 



2n+i> 



(B5) 



where c n {u) are the same coefficients as used in definition 
of M v (z) and T v (z). One can show that W u {z) differs 
from M v (z) only by some constant prefactor. We uti- 
lize this fact calculating the prefactor at some moderate 
values of the argument (z ~ f ), and then evaluate m us- 
ing the function W v (z), that has well defined behavior at 
z -> 0. 

Finally the radial wave function 'J; (r) for partial wave 
I can be written as linear combination of M v (z) and T v (z) 



+ cos(0 ; + § v+ § ) 



T„ In 



, (B6) 



where 0; is a parameter which appear in the small r ex- 
pansion ©. Now, the behavior of ^(r) at small and 



APPENDIX C: EXPANSIONS FOR SMALL 
ENERGIES 



Behavior of the parameters m and v at small energies 
is given by the following series expansions 



4 (i-i)(| + |)(i+3) 



+ 0(<7 4 ), (CI) 



and 



m(q) 



q^l+i r(|-z) 

4 







+ o( g /+5 / 2 ) 



(C2) 



where VK 3 -) denotes the digamma function [5(J |. To obtain 
the above formulas, one can in fact solve the recurrence 
relation (IA4j) keeping only the lowest-order terms Co, c_i 
and c\. 
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